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Abstract. We investigate theoretically the performance of two-qubit resonant gates 
in the crossover from the strong to the ultrastrong coupling (USC) regime of light- 
matter interaction in circuit QED. Two controllcd-PHASE (CPHASE) gate schemes — 
which work well within the rotating wave-approximation (RWA) — are analyzed while 
taking into account the effects of counter-rotating terms appearing in the Hamiltonian. 
Our numerical results show that the fidelity of the gate operation is above 96% when 
the ratio between the coupling strength and the resonator frequency, g/u> r , is of about 
10%. Novel schemes are required in order to implement ultrafast quantum gates when 
increasing the ratio g/oo r . 
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1. Introduction 

Circuit quantum electrodynamics (QED) [TJ [21 E] has become arguably a prominent 
solid-state architecture for quantum information processing [I]. This technology is built 
from superconducting qubits that interact with microwave fields in one-dimensional 
resonators. Its potential applications include single- and two-qubit quantum gates 
via dispersive [5j El [TJ [8] or resonant coupling [91 [101 [HI [12], three-qubit gate and 
entanglement generation [131 EH E3 EEl EH EE] , as well as novel quantum phenomena 
in circuit QED such as dynamical Casimir effect [19], and the ultrastrong coupling 
(USC) regime of light-matter interaction [201 1211 [221 1231 [21]- None of them is accessible 
in existing cavity QED technology. In particular, the prospects of reaching the USC 
regime, where the coupling strength g becomes an appreciable fraction of the resonator 
frequency uj r (0.1 <g/uj r < 1), brings us the possibility of speeding up quantum gate 
operations at time scales of the order of a few nanoseconds. This will help to overcome 
decoherence and implement scalable quantum computation. Nonetheless, this comes 
at the expense of losing the intuitive dynamics of the Jaynes-Cummings (JCM) model. 
Indeed in the USC regime, the well-known rotating-wave approximation (RWA) fails 
as the counter-rotating terms become relevant [25]. This means that in order to make 
possible fast high-fidelity quantum gates, a revision of schemes proposed within the 
RWA is in order. 

In this work, we study two schemes for resonant CPHASE gates in the USC regime, 
which are designed assuming that the RWA holds. In Sec. |2j we analyze the CPHASE 
gate scheme proposed in Ref. [10] considering ratios g/co r in the USC regime. In Sec. [3j 
we propose a modification of the previous CPHASE gate protocol and investigate how 
it performs when increasing the ratio g/u r . Finally, in Sec. |1J we discuss the physical 
platform where these protocols can be implemented, and we present our concluding 
remarks. 



2. Resonant CPHASE gate: scheme I 



In the computational basis of two qubits {|<7i,<72), \gi,e2), \ei,g2), ^1,^2)}, a general 
CPHASE gate is described by the unitary transformation 
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where |^) and | e^) {i = 1,2) are the ground and excited states of the i-th qubit. In 
particular, if 9 = tc, the above unitary transformation leads to a controlled 7r-phase 
gate, that together with single-qubit rotations form a set of universal gates for quantum 
computation [I]. 

In this Section, we study the performance of the scheme proposed in Ref. [10] for 
different values of the interaction strength between qubits and resonator field. This 
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Table 1. Operation steps of the protocol I. 



Step 


Transition 


Coupling 


Pulse 


(i) Mapping 




e 2 ,0) ->■ -i 


02,1) 


9 92,^2 


tt/2 


(ii) CPHASE 




|ei,l) -> - 


ei,l) 


9e lt ai 


7T 


(iii) Back Mapping 




-> -i 


e 2 ,0) 


9g2,e2 


3tt/2 



protocol is based on the resonant interaction of three-level superconducting qubits and 
a single mode of the resonator. In this setup, we assume that the qubit transition 
frequency can be tuned in order to switch selectively on and off its coupling to the 
resonator. The logical qubits are encoded into the two lowest energy levels , |ej), 
while the third state |dj) is used as an auxiliary level. 

This first protocol is displayed schematically in Table [T] and Fig. [1} Let us suppose 
that the initial state of the system reads 

|*in) = (h\gug 2 ) + b 2 \g h e 2 ) + 6 3 |ei,</ 2 ) + &4|e 1 ,e 2 )) ® |0), (2) 

where 6j (i — 1, . . . , 4) are the arbitrary complex coefficients. Within the realm of cavity 
or circuit QED, it is often the case that the interaction between qubits and resonator 
field can be described using the JCM, where the coupling strength is small enough so 
that the RWA is applicable. In such the state evolves to 

l^out) = (b 1 \g 1 ,g 2 ) + b 2 \gi,e 2 ) + h\e 1 ,g 2 ) - h\e 1 ,e 2 )) <E> |0), (3) 

leading to a controlled 7r-gate operation. 

A natural question is whether this protocol can be extended to higher values of 
coupling strength, where the interaction lies in the USC regime. This is a relevant 
question for quantum computation, as it would lead to faster gate operations. We have 
analyzed the fidelity of the above protocol considering the USC regime, where the RWA 




Figure 1. Schematic of the protocol I for a resonant CPHASE gate. 
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breaks down. In this case, The Hamiltonian of the system reads 

i=l,2 

#NRWA = Yl E * \ti ^ I + hbJr 
j=9,e,a 

+ % ei ,ai<,a a (a + a f ) + hg g2)e2 a x g2 ^ e2 (a + a f ), (4) 

where o£ l = \l)(k\ + \k)(l\, Ej i is the energy of the j-th level for the z-th qubit, and the 
9ei,ai, 9g2,e 2 are the corresponding coupling strengths. 

Taking the state ^ as the ideal one to compare with the resulting state of the 
protocol when including counter-rotating terms in the qubit-resonator interaction, the 
fidelity can be computed as 

F — I (^rwa I ^nrwa) 1 2 

= I (*in I ^RWA^NRWA I *in) f \ (5) 

where J/rwa, I^rwa) and C/nrwA) I^nrwa) are the evolution operator and final state in 
RWA and non-RWA cases, respectively. 

For the sake of simplicity, we assume the coupling strengths between each qubit and 
the resonator are equal, g ei>ai = g 92 ,e 2 = g, and study the case of an initial maximally 
entangled state 

l*(o) = ^ (\9u92) + \9i,e 2 ) + \ei,g 2 ) + |ei,e 2 )) ® |0). (6) 

The fidelity of the operation as a function of the normalized coupling strength g/u r is 
shown in Fig. [2] Our simulation shows that the fidelity of the gate operation decreases 
while increasing the coupling strength. For a ratio g/u) r = 0.065, the fidelity drops 
below 0.99, while for g/u r = 0.12 — which was reached in recent experiments [221 123] 
the fidelity goes down to F ~ 0.968, and for g/u r = 0.2, the fidelity is only F ~ 0.89. 
These results mean that, although this protocol could be still used with state-of-the-art 
circuit QED technology (221 E3], its fidelity drops gently as the coupling strength is 
increased beyond these values. 



3. Resonant CPHASE gate: scheme II 

In this Section, we propose an alternative scheme based on a similar configuration, 
but adding an external microwave field driving the transition \e\) -H- |ai) of the first 
qubit, which will assist the gate operation besides the qubit-resonator interaction. The 
advantage of this modified protocol is that it does not require the accurate adjustment 
of the qubit transition frequency. The different steps of this scheme are depicted in 
Table [2] and Fig. [3] — where a single subsystem has a three-level structure. 
Within the RWA, protocol II produces the finial state 

\^ out) = (b 1 \g 1 ,g 2 ) +b 2 \gi,e 2 ) -b 3 \e 1 ,g 2 ) + h\e 1 ,e 2 )^) ® |0), (7) 

that has a 7r-phase shift on state \ei,g 2 ). To go beyond the RWA in the scheme, we must 
take into account counter-rotating terms of both qubit-resonator interaction as well as 
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Figure 2. Fidelity versus normalized coupling strength gju) r for the CPHASE gate 
in scheme I with maximally entangled initial state. 



Table 2. Operation steps of Protocol II. 



Step 


Transition 


Coupling 


Pulse 


(i) Mapping 


|e 2 ,0) -> -i \g 2 , 1) 




tt/2 


(ii) Rotate qubit 1 


ex) — > —i a-i) 


^ei,ai 


tt/2 


(iii) CPHASE 


1) ->■ - l^i, 1) 


9gi,e\ 


7T 


(iv) Back Rotate qubit 1 


ai) — > —i ei) 


^ei,ai 


rr/2 


(v) Back Mapping 


\g 2} 1) -)■ -i |e 2 ,0) 


9g2,e2 


tt/2 




Figure 3. Sketch of protocol II for a resonant CPHASE gate. 
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Figure 4. Fidelity versus normalized coupling strength gju) r for the CPHASE gate 
in scheme II with maximally entangled initial state. 



semiclassical model to describe the qubit and classical microwave field interaction. The 
Hamiltonian for semiclassical model — in the Schrodinger picture — without RWA reads 

#nrwa = E ei \e 1 )(e 1 \ + E ai \ai){ai\ + Ml ei>ai <, ai (e iw * + 4 ), (8) 

with £l ei>ai = d euai £o/fo being the Rabi oscillation frequency, and ojl, £o being the 
frequency and amplitude of the classical driving field respectively. Likewise, the 
corresponding Hamiltonian in rotating frame with respect to qubit frequency one reads 

#nrwa = ^+ (e- 4 ^-^-!)* + e^-m-x-i)*) (9) 

+ Vl e a <r_ (e'^ Wi-Wei ' 01 ^ + e~ i ( UJL+UJei ' ai ^ t ) 

where uj ei>ai = E ai — E ei is the transition frequency between levels |ei) and \ax). 
The wave function in rotating frame can be written as 

|*' (t) ) =C ei (t)\e 1 ) + C ai (t)\a 1 ), (10) 

where C ei (t) and C ai (t) are the complex coefficients for the excited state and auxiliary 
state, respectively. The corresponding evolution equation for these amplitudes are 

C ei (t) = -^^(e^^-^-i^ + e-'^^^)*)^^) (11) 

C ai (t) = - ifi ei ,a 1 (e~ i(&,i ~ Wei,ai) * + j^"*^) C ei (t). (12) 

We have analyzed numerically the protocol fidelity as a function of the ratio g/u r . 
As shown in Fig. |4j it is clear that the fidelity decays in a similar fashion as the previous 
one, as the normalized coupling strength increases, assuming that ^ei,ai = 9- Thus this 
protocol is also unsuitable considering coupling strengths of about g/u r > 0.15, well 
within the USC regime. 

4. Discussion and Conclusion 

Nowadays, the controllability of superconducting devices made of Josephson junctions 
has led to the implementation of different types of superconducting qubits described by 
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phase, charge, or flux degrees of freedom |26j. This fast growing technology now allows 
to access the coherent control of these artificial atoms in a three-level configuration. 
Leading experiments have shown this ability with transmons [27] , the phase qubit [12] , 
and the flux qubit [28] . These advances support the implementation of our protocols in 
various kinds of superconducting circuits. 

Besides, the complexity of the protocols is manifested by the degree of difficulty 
in manipulation, i.e. the number of the pulse sequences needed to perform. Although 
in our protocols three or four steps are needed to perform the gate, and an external 
field is used to assist the operation, we do not require direct qubit-qubit interaction 
as compared with the two-qubit algorithms demonstrated in Ref. [8]. In that scheme, 
the qubit-qubit coupling is realized by a virtual excitation of an intracavity field and 
this second-order based coupling results in a slower operation. In another experimental 
realization [12], the two-qubit controlled-Z and controlled-NOT gates were realized by 
making use of an extra capacitor mediating the coupling between qubits. 

In conclusion, we have analyzed numerically the performance of those two CPHASE 
gate schemes in USC regime, which are designed to work ideally in the conventional 
strong coupling cases. Our results indicate that these two schemes can still work with 
high fidelity for values of normalized coupling strength of about g/w r = 0.12. This is 
the maximum value of coupling ever achieved between quantum light and matter, and 
was realized in recent experiments using superconducting flux qubits j22j [23]. However, 
our numerical analysis also shows that given the present prospects to go beyond this 
coupling strength using similar technology [251 [291 E0] , there is a need to develop new 
protocols for quantum gates beyond the RWA, making possible the design of ultrafast 
quantum gate operations for quantum information processing [3T] . 
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